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Abstract 

We present a new method for reconstructing the cosmological den- 
sity, pecuhar velocity and peculiar gravitational potential on large 
scales from redshift data. We remove the distorting effects of line- 
of-sight peculiar motions by using the linear theory of gravitational 
instability, in which the potential is the solution to a linear partial 
differential equation first derived by Kaiser. We solve this equation 
by deriving its Green function; the Green functions for the peculiar 
velocity and the real-space density follow directly. Reconstruction of 
the cosmic fields is thus reduced to integration over redshift space with 
the appropriate Green function kernel. Our algorithm has a single in- 
put parameter, (3 = Vf'-^ /b, where b is the linear bias factor and 
is the cosmological density parameter. We discuss the virtues of this 
method for error control, for estimating /3, and for constraining the 
bias mechanism. 
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1 Introduction 



Rcdshift surveys are testimony to the richness of structure in the universe 
on large scales (see, for instance, Giovanelli & Haynes 1991). The observed 
mass distribution in "redshift space" , the space defined by two sky coordi- 
nates and a radial recession velocity, depends on both the real-space mass 
density and the peculiar velocity field, each of which are of great impor- 
tance in cosmology. The purpose of this paper is to present a new method 
for evaluating both the density and velocity fields, as well as the peculiar 
gravitational potential, directly from redshift data. 

While quantifying the redshift-space distortions is difficult in general. 
Kaiser (1987) has solved the problem for large scales (> lOMpc) in the 
linear regime of gravitational instability theory (Peebles 1993, §5): The 
effect of linear peculiar flows on the observed redshift-space density field is, 
loosely speaking, to enhance the amplitude of fluctuation modes along the 
observer's line of sight by a factor of 

/3 = J^o-V6, (1) 

where the parameter b admits the possibility that there is a linear bias in the 
fluctuations of luminous matter relative to the total mass distribution. This 
dependence on has been exploited to estimate the cosmic mass density 
from redshift data (Hamilton 1992, 1993; Fisher, Scharf & Lahav 1994; Cole, 
Fisher 8z Weinberg 1994; Bromley 1994). For the purposes of removing 
redshift space distortions, we will assume that (3 is known or that it can be 
estimated in some self-consistent way from the data. 

The problem of reconstructing cosmic fields from redshift data has at- 
tracted considerable interest in the last few years. Yahil et al. (1991), Taylor 
& Rowan-Robinson (1993) and Gramann, Cen & Gott (1994) have imple- 
mented iterative algorithms for real-space density and peculiar velocity re- 
constructions in the linear or quasi-linear regimes. In rough terms, the itera- 
tive strategy is to assume that the redshift data approximate the real-space 
density, calculate the peculiar flows, then correct the density accordingly. 
The procedure is repeated until the real-space distribution and the velocity 
field converge to a solution that is consistent with the observed redshifts. A 
different strategy has been to expand all fields in spherical harmonics, and 
then estimate each multipole individually (Fisher, Scharf Sz Lahav 1994; 
Nusser &: Davis 1994). 

Here we start with Kaiser's equation for redshift space distortions and 
solve for the Green functions that give the real-space density, the peculiar 
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velocity field, and the peculiar gravitational potential. Thus, the reconstruc- 
tion algorithm simply involves integration over the redshift-space density 
with a Green function kernel. The method is designed to hold on scales 
in the linear regime under the assumption that light traces mass or that 
luminous matter is a linearly biased sample of the total density field. The 
technique is valid for any survey geometry. In particular, we do not make 
Kaiser's "plane-parallel" approximation wherein the redshift survey is as- 
sumed to lie in a compact region at a large distance from the observer (see 
Zaroubi & Hoffman 1994). 

The Green function method has distinct advantages over previous tech- 
niques: it is a noniterative, linear algorithm applicable to any survey geom- 
etry, and the work is done entirely in real space, not in the Fourier domain 
or with spherical harmonic coefficients. Hence the error bars due to edge 
effects can be calculated in a straightforward and intuitive way. 

This paper is organized as follows. After reviewing Kaiser's result and 
establishing some notation in Section 2, we derive the Green function for 
the peculiar gravitational potential in Section 3, from which the kernels for 
peculiar velocity and real-space density follow directly. In Section 4, we spe- 
cialize to the simple case known as the "distant survey approximation" . In 
Section 5, we discuss how the technique can be applied to real-world situa- 
tions where one faces nonlinear clustering, shot noise and finite volume. In 
Section 6, we give a numerical example of density reconstruction from artifi- 
cial redshift-space data. We conclude by summarizing possible applications 
of the method for extracting Q and constraining the biasing mechanism. 

2 The partial differential equation 

The first step is to obtain an expression, analogous to a Poisson equation, 
which relates the real-space peculiar gravitational potential (p to the density 
fiuctuation field 5s in redshift space (we denote redshift-space quantities 
with a subscript-^ while subscript-i? indicates real space). This was first 
done by Kaiser (1987), but we give a brief derivation here for completeness. 
The fluctuation field is defined in the usual way, 

S=[pir)-p]/p, (2) 

where p is the mass density and p is its mean. Here, the density fields are 
associated with luminous matter, and we assume that fluctuations in light 
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are proportional to fluctuations in total mass on the scales of interest. We 
use time units such that the Hubble constant Hq = 1. 

Denoting the radial part of the peculiar velocity field v at a real-space 
position r by 

t/(r)^v(r)-f, (3) 
and the radial peculiar velocity in the observer's frame by 

AC/(r)^[v(r)-v(0)]-f, (4) 

the apparent position s of a galaxy in redshift space is 

s(r) = r + AC/(r)f . (5) 

As long as the perturbations are sufficiently small that this mapping is one- 
to-one, the overdensities in redshift space and in real space are related by 

(1 + 6R)d^r = (1 + 6s)dh . (6) 

Evaluating the Jacobian of the mapping (|5|) yields 

d^s= (^l + ^^\l + drU)d^r, (7) 

where dr denotes the radial derivative d/dr. Thus, to first order in v and 
5, and in the absence of selection effects and survey limits, the equation 
relating the densities in redshift space and in real space is 

6s = 6r- drU - 2- + 2 ^^°^ . (8) 

r r 

We can express the first three terms on the right hand side in terms of the 
peculiar gravitational potential, cj). The density Sr and cj) are related through 
the Poisson equation, 

= A5r, (9) 
where (cf. Peebles 1993, eq. [5.107]) we have introduced the constant 

A = AttGpo = (10) 

In linear theory, the peculiar velocity is simply proportional to the gravita- 
tional force field, and is given by (cf. Peebles 1993, eq. [5.117]) 

v = -^V<A. (11) 
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Thus the velocity terms are 



drU=-^d',<P; J = -^^dr^- (12) 

Substituting this into equation (^) yields what we will refer to as the Kaiser 
equation, 

Del) = AS's , (13) 
where the differential operator D is defined as 

D = V'^ + Pd^ + 2(3-dr . (14) 
r 

Here we have eliminated the last term in equation by defining 5'^ as the 
redshift space density observed in the cosmic rest frame, the frame relative 
to which V is defined. From the observed dipole anisotropy in the microwave 
background radiation, we know that the peculiar velocity of our solar system 
is (Smoot et al. 1992) 

v(0) = 370 ± 10 km/s (15) 

in the direction a = 11.2'', 5 = —7°. The cosmic variance in this dipole is 
merely a few km/s, i.e., considerably smaller than the errors quoted above. 
S'g can be directly computed from the redshift-space density we observe, 6s, 
either by using the first order correction derived above, 

6's{r) = 6s{r)-2^^^, (16) 

or by simply adjusting the observed redshifts (radial galaxy positions): 

s' = s + v(0)-f. (17) 

In what follows, we will drop the prime and let 6s denote the corrected field. 
It should be noted that this distinction is of only marginal importance. The 
last two terms in equation (^), are usually negligible, and have been dropped 
by most authors. We leave them in because they preserve the self-adjointness 
of the operator D. It turns out that the approximation of dropping them 
gives a more complicated solution, not a simpler one. 
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3 The Green functions 



We now solve the Kaiser equation (13). Comparing the differential operator 
D with the expression for in spherical coordinates shows that it is simply 
the Laplacian with the radial part amphfied by a factor [1 + (5). It is self- 
adjoint for ah /3, which means that its eigenfunctions are orthogonal, and 
that expanding the fields in these functions (spherical Bessel functions of 
non- integer order) provides a straightforward way of solving the equation. 
A related formal solution is given by Nusser &; Davis (1994), but for a 
slightly different operator yielding the solution in redshift space, not real 
space. Here we will instead give a more direct solution, in real space, which 
renders orthogonal expansions unnecessary. 

By the linearity of the problem, we can always write down a formal 
solution as 

^{r) = A J G4r,r')6s{r')dV , (18) 

where the Green function (integral kernel) satisfies 

iDG^){r,T')=6Hr' -r) (19) 

and the Dirac delta function is distinguished from the fluctuation fields 6ji 
and 6s by the absence of a subscript. The constant A in Eq. (|l8|) was 
inserted to allow us to think of G^f) as the operator inverse of D. Although 
Gff) is a function of six variables, there are in fact only two quantities upon 
which the dependence is non-trivial. Since the Kaiser equation is invariant 
under rotation and rescaling, it is easy to see that we can write 

G,(r,r') = «%:^ (20) 



for some function g, with the angle cosine /i = r • r'. Let us utilize this 
property by expanding in Legendre polynomials as 

1 oo 

G^r, r') = + l)P,{^^)ae{r/r'). (21) 



The Kaiser equation implies that the radial functions ai{x) satisfy 

6{x-l) 

4 ' 1 + /? J 47r(l + /3) 



X x-^ 
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Requiring ai to be continuous and imposing the boundary conditions that 
ai{x) remain finite both as x ^ and as x — > oo, this equation has the 
solution 



87r(l + [3)ai 

where the exponent is +a£ if x < 1 and —ai if x > 1, with 



Substitution back into equation (^) gives the solution 

GJv,r') = ^—^T (^^) Pdp) (—X\ (25) 

where r< and r> denote the smaller and larger of the numbers r and r', 
respectively. If we let /3 — > 0, the Kaiser equation approaches the Poisson 
equation, {21 + !)/«£ — > 2 and equation (p5[) reduces to the familiar 

G^M = -^£r^y (26) 

We now consider the peculiar velocity field v, which according to equa- 
tion ([TT| ) is proportional to the gradient of the peculiar gravitational poten- 
tial (j). Rather than computing Vi;^ numerically, it is convenient to define a 
Green function Gv expressly for the velocity field which satisfies 



V r 



j Gv(r,r')55(r')dV', (27) 



where, evidently, 

Gv = (5VG^ . (28) 
Taking the gradient of equation ( p5|) with respect to r gives 



where 'je = —1/2 + 0; if r < r' and j£ = — 1/2 — ai if r > r' . The coordinate 
6 refers to a spherical coordinate system with the zenith in the r'-direction, 
i.e., where r • r' = cos 9. 
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For completeness, we finally construct the Green function for the real- 
space density, G^, such that 

5Rir) = I Gsir,r')6sir')d\' . (30) 

Since the Poisson equation gives 

Gs = V^G^, (31) 
taking the Laplacian of equation (p5|) yields 



where the delta- function term represents the full contribution of at r = r'. 

As expected, 

Gs{v,r')^6^{r'-r) (33) 

in the limit /3 ^ 0. Contour plots of Gs and the Green function giving the 
potential, with /? = 1, appear in Figure 1. Note that G(^(r, r') = G^(r',r), 
since is self-adjoint, whereas Gv and Gs lack this symmetry. 



4 The distant survey approximation 

An interesting limit occurs for a distant-survey (Kaiser 1987, Zaroubi &: 
Hoffman 1994) in which r ~ r' or |r' — r| <C r. For definiteness, let us choose 
r' = z, the unit vector in the z- direction. Then dr ^ dz, the radial part of 
the Laplacian ~ d^, and the Kaiser equation reduces to 

dl + d'^ + il + /3)a2] = A5s. (34) 

A simple rescaling of the z-axis transforms the above expression into the 
Poisson equation. This fact and the application of the gradient and Lapla- 
cian operators to the Green function for the potential immediately give us 
the three Green functions in the distant-survey limit: defining 

Ar = r — r' = (sin cos (/), sin sin (j), cos '&)Ar (35) 
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and m = cos(i?), with -& being the angle between Ar and the hne-of-sight 
(the z-axis, parallel to both r and r' in this limit), 

^^^"''') = -4;i^(^ + ^-M~'^'' ^^^^ 
1 / \ -1/2 / \ 

iTTAr^{l + P- Pfify'^ 3 + 3/3 
Thus Gs vanishes on the cone given by 

\z' -z\= + (39) 

takes positive values inside this cone and takes negative values outside of it. 



5 The real world: nonlinearity, shot noise and fi- 
nite volume 

The above treatment applied to an idealized world were the Kaiser equation 
was strictly valid, and where we could measure the redshift space density 
field Ss accurately throughout all space. Alas, this is not the world in which 
we live. Below, we discuss the following three sources of error: 

1. The Kaiser equation, which is exact to first order, is only accurate when 
6 <^ 1 and thus breaks down on small scales where nonlinear evolution 
has become important. 

2. Our galaxy surveys do not measure the continuous field Sr but the 
discrete galaxy distribution, so our estimates of Ss are contaminated 
by Poissonian "shot noise". 

3. Our galaxy surveys sample only a finite volume, which forces us to 
truncate our volume integrals and accept truncation errors in our re- 
constructed fields. 

There are of course many more potential sources of error, such as nonlinear 
biasing and flux errors to mention two, but we will restrict our discussion 
to the three listed problems, which are always present. In what follows, we 
also limit ourselves to reconstruction of the density field, the analogy for the 
velocity field and the peculiar gravitational potential being obvious. 
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5.1 Nonlinear clustering 



On small scales, the gravitational evolution has entered into the deeply non- 
linear regime, where the "fingers-of-God effect" causes virialized clusters to 
appear elongated along the line of sight. Although some authors have at- 
tempted to extend the Kaiser equation into the quasi-linear regime using the 
Zel'dovich approximation, an exact reconstruction of the real-space density 
is of course impossible in the deeply non-linear regime. Essentially, this is 
because the gravitational n-body problem exhibits chaotic behavior, render- 
ing the inversion problem numerically ill-posed. We thus follow the same 
approach as all previous authors, and smooth all our fields {5s, 6r, etc.): 

6s{r) ^ {W*6s){r) = J W{r' - r)6s{r')(fr' , (40) 

the star denoting convolution. For definiteness, let us choose a Gaussian 
smoothing kernel 

and choose the smoothing scale to be large enough that the resulting fields 
can be expected to obey linear theory]^. 



5.2 Shot noise 

What we measure is of course not the smooth field 6s but the discrete galaxy 
density n(r) = J2 ^{^ ~ ^i): with a delta function at the location of each 
galaxy. We estimate the unsmoothed density by simply 6n/n = {n — n)/n, 
where n, the selection function, specifies the survey geometry by giving the 
expected density of observed galaxies at each position in space. Our estimate 
of of the redshift space density 5s is thus the convolution of n/n — 1 with the 
smoothing kernel W. Hence our real- world generalization of equation (]30| ) 
becomes 

" (rO, (42) 



^-1 
n 



where the prime on the left hand side has been added to distinguish the 
estimate 5^ from the true value 5r. We will return to the issue of how to 
choose V , the volume to integrate over, in the following subsection. Since 

^ Since smoothing suppresses high-frequency Fourier components, its success of course 
hinges on the standard assumption that the gravitationaUy induced mode couphng does 
not cause a significant "blue leak" of power from short to long wavelengths. 
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convolution is associative (simply corresponding to multiplication in Fourier 
space), equation (p2|) is equivalent to 



6'^{r)= I {W^Gs){v,v') 
Jv 



n{r') 



(43) 



i.e., to smoothing the Green function (with respect to r') instead of the 
galaxy distribution. 

Let us, as is customary, model n as a Poisson process whose intensity (av- 
erage point density) is A(r) = n(r)[l + (55(r)]. We model (the smoothed) 5s 
as a random field with {6s{r)) = and define its correlation function as 

Ur,r')^{Ssir)6sir')), (44) 

but make no assumptions of it being homogeneous, isotropic or Gaussian. 
Defining the reconstruction error as 

A(r) ^ 5M - 6n{r), (45) 

we now proceed to compute the statistical properties of this error. Since a 
Poisson process satisfies 

(n(r)) = A(r), (46) 
(n(r)n(r')) = A(r)A(r') + 5(r - r')A(r), (47) 

a straightforward calculation shows that 

(A(r)) = 0, (48) 
(A(r)2) = a„(r)2 + at(r)2, (49) 



where 



n(r') 

at{rf ^ = f f Gs{r,r')Gs{r,r")^s{r'y)d^r'd^r". (51) 

Jyjy 

When computing cj^, the shot noise error, the integral is to be taken over 
the volume V (see below). When computing cj|, the truncation error, the 
integrals are to be taken over all space except the volume V. Equation (|48| ) 
simply tells us that our reconstruction is unbiased, in the sense that the 
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expectation value of the error is zero. The remainder of this subsection is 
devoted to the shot noise term in Equation (^9|). 

Although the shot noise error (T„ is easily computed numerically using 
equation (|50|), it is instructive to make a crude estimate of it. First, let us 
assume that n(r) is constant, so that we can factor it out of the integral 
in equation ([50|) . Second, let us make the distant survey approximation for 
G5 (this is a good approximation if r is more than a few smoothing lengths 
away from us, as it is easy to show that the smoothed Green function W-kG 
falls off rapidly, at least as fast as r~^, for |r — r'| ^ r^, so that the main 
contribution to the integral comes from within a few smoothing lengths of 
r). Third, let us overestimate slightly by extending the integration to 
all space (this makes only a minor difference, as we will find it optimal to 
integrate considerably beyond the smoothing scale anyway). Now we can 
apply Parseval's theorem and the convolution theorem to the remaining 
integral, and obtain 

^ (2^ / \Gsm'\^m'd% (52) 

where hats denote Fourier transforms. From equation (|3^), we immediately 
recover Kaiser's familiar result 



Gs(k) 



(53) 



so substituting this into equation (|5^, we obtain the estimate 

1 /• A;2 sine ^, ^^^^ /arctan^\/l\ 

Defining the effective smoothing volume 

Vs ^ (^^^] rl (55) 
Y arctan ^/P J ' 

we can conclude this section with the familiar-looking shot noise formula 

"""'i^ (56) 

In other words, the crucial quantity is the expected number of galaxies in a 
smoothing volume. 
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5.3 Finite volume 



Although the truncation error at is straightforward to compute numerically 
given any prescribed volume V, either by evaluating the integral in equa- 
tion ( ^T|) or by making reconstructions from Monte Carlo galaxy catalogs, 
we now make a very crude estimate at to obtain a qualitative understanding 
of how to best choose V and the smoothing scale. Although the unsmoothed 
fields 6r and 6s have ample structure on very small scales, the correspond- 
ing smoothed fields are virtually featureless on scales r <C r^. Hence the 
correlation is almost perfect for |r' — r| <C r^, and typically falls off like 
some power law for |r' — r| ^ Vg- Thus let us make the crude approximation 

for some constant a which is the r.m.s. fluctuation in a smoothing volume. 
The only non-negligible contributions to the integrals in equation ( [5l| ) thus 
arise when |r' — r| < r^. This can be interpreted as there being a large 
number of independent volumes in which the field is coherent, and that the 
total variance at is simply a sum of the variance arising from each coherence 
volume that we do not look at, i.e., that falls outside of V. Assuming that 
V is large relative to the smoothing scale, replacing ^^(r, r") by ^^(r, r') 
will thus be a good approximation, and equation (|5ll) reduces to simply 

at{rf « aa^Vs f Gs{v,v'fd\\ (58) 

where a is some constant of order unity. In the distant survey approximation 
of equation (38), we see that Gs{t, r') oc |r' — r|~^ far away, so if V is chosen 



to be a sphere of radius Vc around r, then 

5.4 How to choose V and the smoothing scale 

In this section, we give some useful rules of thumb for how to choose the 
volume V and the smoothing scale Vg. 

Since the integrand in equation (|5^) is always positive, the shot noise er- 
ror an always increases if we increase the volume V , whereas equation ( |59| ) 
shows that the truncation error at decreases as we increase V . Since all 



12 



we care about is the total error cj^ + af, we find ourselves facing a trade- 
off. Similarly, equation (|5^) shows that the shot noise error decreases if we 
increase the smoothing volume Vs, whereas equation ( ]59| ) shows that the 
relative truncation error ctj/ct increases if we increase the smoothing volume 
Vs- Minimizing the total error by differentiating cr^ + with respect to 
V and Vs, we of course end up with parameter choices such that the two 
sources of error are of comparable importance, i.e., such that (T„ is of the 
same order of magnitude as at- Another way of phrasing this is that given 
some finite level of truncation error at, there is no point in trying to re- 
duce the shot noise error an way below this value (and vice versa), as this 
will cause almost no reduction in the total error. We thus obtain one of 
our rules of thumb by simply equating an and at, which after substituting 
equations ( ^6|) and (|59|) gives 




so that the optimal smoothing volume Vs then goes roughly as the geometric 
mean of the total volume V and the average volume per galaxy, 1/n. 

In practice, the survey specifics often place a firm lower limit on either 
cr„ or at, which makes the best choices of V and Vg quite simple. We now 
turn to a couple of specific examples: 

5.4.1 Slices and pencil beams 

If the survey geometry is that of a thin slice, such as for instance the CfA or 
Las Campanas redshift surveys, or that of a so-called pencil beam, then the 
dominant source of error at a point r will usually be the truncation error 
at caused by the most nearby edge of the slice/pencil. This means that 
we should choose the smoothing scale as small as we feel comfortable with 
given the non-linearities on small scales, perhaps around ~ 10/i~^Mpc. 
Equation (^8|) shows that since Gs'^ falls off so rapidly (about as |r' — r|~^), 
the truncation error will be completely dominated by the smallest distance 
from r to the edge of the volume V. In other words, given any volume V, 
we can get almost as small truncation errors by replacing it by the largest 
sphere around r that will fit inside it, and this defines a characteristic at 
that we refer to below. In practice, the reconstruction calculations are so 
fast that there is little need for such attempts to save computer time, espe- 
cially considering the amount of effort that has gone into the preceding data 
collection. Rather, the main concern is that we want to avoid extending 
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V out to such large radii that the shot noise error becomes as large as at- 
Since n(r) is usually a decreasing function of r, Equation ( |56|) shows that 
this would happen if r where allowed to be so large that n(r) < l/atVg- In 
conclusion, when reconstructing the real-space fields at a point r, we will do 
fine if we simply choose V to be the entire slice volume, but truncated at 
some cutoff radius Vc such that 

• rc—r exceeds the sideways distance from r to the edge of the slice/pencil 
beam 



5.4.2 "Non-skinny" surveys 

For surveys where all three dimensions are large relative to the smoothing 
length (such as the IRAS 1.2Jy all-sky survey), the situation is quite dif- 
ferent. When reconstructing the real-space field at points r that are many 
smoothing lengths away from the survey boundaries, truncation errors will 
not a priori dominate as above. Hence an efficient analysis strategy is to 
simply truncate the data at some radius Tc where n begins to fall off rapidly, 
take V to be the the entire remaining survey volume (since an saturates a 
few smoothing lengths away from r anyway), and then select the smoothing 



scale according to equation (60) 



6 Examples 

We now evaluate the Green function method in two simple numerical exam- 
ples. First, for a qualitative assessment, we realize a single 3D field in real 
space, remap the density into redshift space with peculiar fiows derived from 
linear theory, and then perform the reconstruction of the real density field 
with the Green function method. The results of this procedure are given in 
Figure 2 for an Q = /3 = 1 Universe, where 2D slices of the density fields are 
shown. In addition to the reconstructed field (lower-left plot) we provide 
a reconstruction using an incorrect /? value of 0.25 (lower-right plot). This 
illustrates the sensitivity of the method to /3, a feature which may be useful 
in constraining and the degree of linear bias. 

Details of the procedure that led to the data in Figure 2 are as follows: 
We start with a triply periodic 250 Mpc cube and realizes a Gaussian 
random field on grid of 64^ points. The power spectrum of density fiuctu- 
ations has the form of standard cold dark matter (e.g., Efstathiou, Bond 
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& White 1992), smoothed with a Gaussian filter on scales of 50 h"^ Mpc. 
The normalization is chosen near the COBE value to approximate observed 
large-scale structure; before smoothing, the r.m.s. density fluctuations in a 
sphere of radius 8 Mpc are at a level of 1.1 relative to the mean density, 
i.e., (Tg = l-l- We then estimate the velocity field at the grid points and 
displace the the points by an amount appropriate for an observer at some 
chosen origin to mimic the distribution in redshift space. These points are 
interpolated onto a lower resolution grid (32^) to reduce discreteness effects. 
This lower-resolution grid represents the volume-limited red-shift space field 
that is used in the reconstruction. For the 2D plots in Figure 2, we perform 
the reconstruction on a plane of grid points; at each point a 3D integral 
over the redshift-space grid is evaluated with the kernel of equation p^ ) 
under the assumption that the density field is a set of weighted Dirac-delta 
functions on the grid. For speed enhancement, the sum in equation ( |32|) was 
evaluated at a grid of points (x, fi) once and for all in a lookup table. When 
r ~ r', the order zero analytic reconstruction kernel (Ozark) of equation ( |38| ) 
was used. 

Evidently, the reconstruction works well for regions away from the edges 
of the mock survey in Figure 2. The reconstruction has recovered both the 
amplitude of the original field and the topology of the isodensity contours. 
However, a qualitative assessment of the errors is in order. Our focus in this 
section is on truncation errors as well as errors incurred during interpola- 
tion and integration. We estimate the errors with a Monte Carlo method; 
we realize 72 random density fields, perform reconstruction on each as de- 
scribed above, and determine the r.m.s. point-wise difference between the 
real density field and the reconstruction. Figure 3 shows the r.m.s. error in 
reconstruction for points on one of the principle axes of the cubical survey 
volume (the points lie on a path which bisects the plane shown in Figure 2). 

Truncation effects are clear in Figure 3 from the increase in error near 
the boundary of the survey volume. The truncation error worsens as the 
correlation length of the density field increases, as can be seen from a com- 
parison of the error curves with correlation lengths of 50 and 100 Mpc 
(induced by smoothing on those scales) ; large correlation lengths correspond 
to low frequency modes which are more susceptible to the ill effects of a finite 
survey volume. 

For the case of the realizations with a smoothing scale of 50 Mpc, 
we estimate the theoretical truncation error by directly integrating equa- 
tion (|5T|). The relatively high dimensionality of the integral and the rapid 
variation of the integrand near the boundaries of the survey volume suggest 
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that this problem is best handled by parallel computation (we are grateful 
for time on a Cray T-3D at Jet Propulsion Laboratory). For reference, an 
approximate theoretical error curve is shown in Figure 3. The remaining er- 
ror comes from grid interpolation, the numerical integrator, and to a lesser 
extent the breakdown in the Kaiser equation, which is of course only exact 
for (5 ^ 1. This latter source of error scales roughly as (5^, and does be- 
come more important at smaller smoothing scales. While in this example 
we are content to let these errors lie between 5% and 10% near the center 
of the survey, they may be further reduced with a finer grid mesh and a 
higher-order integrator. 

7 Discussion 

We have presented the Green functions for computing the gravitational po- 
tential, the peculiar velocity field and the real-space density fiuctuations 
from the distribution of matter in redshift space. Our results are based on 
Kaiser's (1987) analysis of redshift-space distortions in the linear regime, 
and are applicable to surveys of arbitrary geometry. 

One virtue of the method presented here is that the entire calculation is 
performed in real space. This means that the effects of complications such as 
survey boundaries and the presence of galaxy clusters, which are manifestly 
in physical space rather than Fourier space, are easy to understand and 
quantify. For example, the magnitude of the Green function in the parts 
of space that we are ignoring immediately gives an indication of how much 
information we are loosing, i.e., on the magnitude of the error bars on the 
reconstructed field. As we have discussed, this allows the problems caused 
by shot noise, survey boundaries and nonlinear clustering to be minimized 
by appropriate choices of the smoothing scale and the integration volume 
V. 

We wish to emphasize several possible uses of the Green function method. 
The first is toward ferreting out estimates of the linear growth parameter 
/? = OP-^/b. While /3 has been estimated in elegant manners from the cor- 
relation anisotropy (e.g., Hamilton 1992), it may also be determined with 
the Green function Gs- The best estimate for (5 is the value which mini- 
mizes the anisotropy in the real-space two-point correlation function of the 
reconstructed field. Quantitative conditions for minimizing the anisotropy 
in may be derived with an orthogonal decomposition in spherical har- 
monics (e.g., Hamilton 1992). For example, a least-squares criterion may be 



16 



imposed on the amplitudes of terms with spherical harmonics Yim for which 
I > 0, with greatest weights assigned to quadrupole terms, as they couple 
strongly to the redshift-spacc distortions. 

The second application of the Green function method is to extract infor- 
mation on the nature of biasing. Specifically, the reconstructed real-space 
density field contains the implicit assumption of a linear bias between lu- 
minous and dark matter. Thus it is worth comparing the field with den- 
sity reconstructions from other techniques such as the POTENT algorithm 
(Bertschingcr &: Dekel 1989). While such a comparison will not alone break 
the degeneracy between b and J7, it can test the validity of the linear biasing 
hypothesis. 

Ultimately, the method's most important use is the estimation of peculiar 
velocities. Even linear peculiar flows can provide stringent constraints on 
cosmological models (e.g., Lauer &; Postman 1994; Tegmark, Bunn Sz. Hu 
1994). Since the Green function method generates peculiar velocities directly 
from redshifts without the difficult measurement of real-space distances, it 
may have fruitful application in the large-scale digital sky surveys that are 
currently coming on line. 

Our purpose here has been to provide the mathematical underpinnings 
of the Green function approach to cosmic field reconstruction from redshift 
data. We are now embarking on the next step, the analysis of redshift 
surveys. 

We thank James Anglin, Marc Davis, Andrew Hamilton, Saleem Zaroubi 
and our referee for useful comments and suggestions. BCB acknowledges 
support from the NASA HPCC program and help from Isabel Dulfano. 
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Figure 1: The Green functions and Gs, for (3 = 1. 
To the left of the y-axis are contours of G^(r, r'), while the contours to 
the right correspond to Gs{r,r') for r = {x,y,0) and r' = (0,1,0). The 
light solid lines and dashed contours indicate positive and negative values, 
respectively, with logarithmic spacing. The heavy solid line is the Gs{to) = 
contour. The trend in both functions is toward high absolute values near the 
center of the figure. Values of these functions for more general choices of r 
and r' follow directly from symmetry and scaling properties. The functions 
essentially give the contribution of a point mass at r to the reconstructed 
field at r' (point at center of the figure) . 



19 



Real Space Density 



Redshift Space Density 




Reconstruction (/S= 1) Reconstruction (/?= 0.25) 



Figure 2: An example of the reconstruction method. 
As discussed in the text, the reconstruction method is applied to a linearly 
evolved Gaussian random field in a 250 Mpc cube. The linearly spaced 
contours, denoting isodensity surfaces at a mid-plane of the cube, show over- 
densities (solid lines), underdensities (broken lines), and the mean density 
(heavy solid lines). The observer is located at a corner of the cube. The 
upper-left plot is the density in real space, generated from a CDM power 
spectrum with 50 Mpc smoothing. The corresponding redshift-space 
density (upper-right) is used to reconstruct the real-space fields; the lower- 
left plot was derived with the correct P value while an incorrect value yielded 
the plot on the lower right. 
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Figure 3: Density reconstruction errors. 
The r.m.s. reconstruction errors, a{r)/^/^^^^{0), are estimated from 72 real- 
izations of a linearly evolved CDM field, with 100 Mpc smoothing (top 
curve) and 50 Mpc smoothing (middle curve), in the survey region used 
for Figure 2. The horizontal axis is distance along a principle axis of the 
cube. The lower (dashed) curve is an estimate of the expected truncation 
noise, (7t(r) for the case of 50 Mpc smoothing. 
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